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Abstract. We study the dynamics of Bose-Einstein condensates in time- 
dependent microtraps for the purpose of understanding the influence of the mean 
field interaction on the performance of interferometers. We identify conditions 
where the non-linearity due to atom interactions increases the sensitivity of 
interferometers to a phase shift. This feature is connected with the adiabatic 
generation of a dark soliton. We analyze the robustness of this phenomenon with 
respect to thermal fluctuations, due to excited near fields in an electromagnetic 
surface trap. 



PACS numbers: 03.75.Lm, 03.75.Dg, 03.75.Gg 



There is currently a large interest in employing Bose-Einstein condensates (BECs) 
of dilute atomic gases in the field of matter wave interferometry [TJ [2J. The 
advantages of atom interferometers compared to optical ones are well-known [TJ [3] : 
greater precision due to the large atomic mass, sensitivity to vibrations, inertial, 
and gravitational forces, access to quantum decoherence and to atomic scattering 
properties, to quote a few. A promising route towards compact applications 
are integrated devices, using for example miniaturized optical elements built near 
nanostructured surfaces [H [5] . 

Atom- atom interactions that play a crucial role in BECs, are usually considered 
as a drawback for matter wave interferometers because they introduce additional 
phase shifts and, more fundamentally, quantum fluctuations and diffusion of the 
relative phase between the parts of a spatially separated BEC, see, e.g., [6]-[13]. In 
this paper, we show that an operation mode for a BEC interferometer exists where 
atom interactions actually enhance the phase sensitivity. We study in particular 
the temporal scheme proposed in Refs. [H [14] where a trapped atom sample is 
split and recombined by slowly deforming the trapping potential. A similar setting 
also describes approximately the flow of a condensate through an interferometer 
with spatially split arms [T5]-[T7]. Interference is usually looked for in the lowest 
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vibrational modes of the recombined trap, whose relative populations depend on 
a phase difference imprinted, e.g., during the split phase (see, however, [T7] for a 
multi-mode interferometer with non- interacting atoms). We show here that when a 
condensate is recombined, a grey soliton is formed and starts to oscillate with an 
amplitude controlled by the interferometer operation. It is quite interesting that 
conversely, the interference phase can also be read out by measuring the amplitude of 
either the soliton oscillation or the condensate dipole mode. In both cases, one achieves 
a better phase sensitivity with detection schemes that are relatively simple compared 
to projective measurements in the vibrational mode basis. Finally, the enhanced phase 
sensitivity of the condensate interferometer is illustrated by numerical simulations 
where the output signal is scrambled due to a fluctuating, random potential. This 
models thermally excited magnetic fields that occur in surface-mounted microtraps, 
and our results allow to identify the noise level that can at most be tolerated. 

In our model for the BEC interferometer, we focus on the quasi one-dimensional 
regime typical for elongated traps, assume zero temperature and describe the 
condensate dynamics along the loosely bound axis by a one-dimensional order 
parameter <&(x,i) that solves the nonlinear Schrodinger or Gross-Pitaevskii equation 



The potential V(x,t) in Eq.([T]) is harmonic for t < and t > r. Similar to the 
setup proposed in Refs.[4lfl4]. a barrier is adiabatically and symmetrically raised and 
lowered during the operation time r, splitting and recombining the condensate. A 
relative phase 6 is imprinted at time r/2, when the two wells are sufficiently separated 
(negligible tunnelling). In the following, we use harmonic oscillator units defined by the 

1/2 

initial trap frequency and adopt ui = [h/(mil)} 1 as unit of length. We normalize 
the order parameter to unity and get a dimensionless effective interaction constant 
g = 2N(a s /ui)(fl±/n) (1 — 1.4603a s /uj^) _1 [19] . where N is the number of atoms, 
a s is the three-dimensional scattering length, Q± is the frequency of the transverse 
(radial) confinement and u± the corresponding ground state size. We assume a 
constant, tight radial confinement and neglect the coupling between radial and axial 
excitations. The numerical integration of Eq.([T]) is done with the split-operator 
algorithm [20j whose results we have checked for convergence using the conservation 
of total energy and the symmetry under time reversal. A typical density distribution 
is shown in FigfTJleft) with a trap potential V(x,t) = \ [x 2 — d(t) 2 ] 2 / (x 2 + d(t) 2 ), 
where the well separation is parametrized by d(t) — 2a sin 2 (nt / t) for < t < r. The 
imprinted phase shift is close to tt so that in the ideal gas case, one expects an output 
state close to the first excited trap eigenstate. One clearly sees that a small deviation 
from this phase shift gives a grey soliton that oscillates in the harmonic trap. 

A more quantitative characterisation is usually based on the final populations po,i 
of the lowest trap eigenstates (k = 0, 1) 



where |$(t)) is the order parameter at time t = r. We choose here the convention 
that |</»o, l) are the ground state and the first excited state in the harmonic trap, 
including the nonlinearity. We illustrate below that this improves the accuracy of a 
two-mode approximation (see Figj2j) . Two alternative quantities, that appear simpler 
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Figure 1. (left) Recombining the interferometer arms after a phase shift close to 
tt produces a grey soliton oscillating in the harmonic well (dark zigzag lines), as 
well as dipole and breathing oscillations of the background condensate (bright). 
Parameters chosen in the numerical simulation: g = 10, G = 0.9 w, operation time 
t = 70, maximum splitting 4a = 8. (right) Evolution of the weights |co,i(t)| 2 
for the first two eigenstates of the GPE in the instantaneous trapping potential 
V(x, t) (top and centre) and of the weights of higher modes (bottom). 



to be measured in practice, are the condensate mean position after recombination and 
the centre of the density dip characteristic for the grey soliton 



{x) t = /dx|$(x,i)| 2 x, 



q t = miniocai|$(:M)| 5 



for t > t. (3) 



These quantities show oscillations whose amplitude depends on the imprinted phase 
shift O, as shown in Fig. |2fleft). The key observation is that around a complete 
swapping to the odd output state, the curves are much narrower compared to the 
non-interacting case where the dipole mode oscillation amplitude is | sin Q\/V2. In 
this sense the nonlinearity increases the phase sensitivity of the interferometer^ 
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Figure 2. (left) Maximum oscillation amplitudes of the mean position of the 
whole condensate and of the soliton position, as a function of the phase shift ©. 
For phase shifts very different from © = tt, the grey soliton reaches the condensate 
border in the recombination stage and escapes; for these values, no results are 
plotted. Mean position amplitude (x6 for clarity): g = (dash-dotted line), 
g = 5 (thin dashed line with triangles), g = 10 (thin solid line with squares). 
Soliton amplitude: g = 5 (thick dashed line with dots), g = 10 (thick solid 
line with stars), (right) Final population po of the trap ground state, computed 
by solving the full nonlinear Schrodinger equation (thin solid line with circles, 
g = 10), the two-mode model (solid line, g = 10) of Eq. (5), and the ideal Bosc 
gas (dashed line, g = 0). 

§ We have checked that the soliton oscillation agrees with the effective equation of motion of Ref . 121 1 
in the Thomas- Fermi limit, i.e. for interactions g ^> 1. The difference with respect to Ref. 1221 may 
be only apparent and related to the different scaling convention for the dimensionless parameters. 
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The populations po,i show a similar behaviour which can be understood using a 
two-mode model similar to those studied in Refs.[T31[T51[IE]- We consider an expansion 
of the order parameter 

${x,t) = Y J Ck{t)<l>k{x\d{t)), (4) 
k 

where cj)f.(x; d(t)) are (real) eigenstates of the GPE in the double- well potential with 
d(t) and g fixed. In the adiabatic approximation and restricting the expansion to the 
lowest even and odd modes k = 0, 1, the equations for the coefficients are 

ic k = [i k c k +#(2£>oi - O k k) \ci- k \ 2 c k + gO m c\_ k c* k (5) 

where fi k are the chemical potentials of the eigenstates. Due to parity conservation, 
the only nonzero interaction matrix elements are 

o kl = Jdx[M^d)Mx;d)} 2 . (6) 

Starting from the ground state 4>o(x; d(0)), adiabatic evolution and phase imprint 
lead to the amplitudes c (t/2 + 0) = cos(9/2)co(r/2 - 0) and ci(r/2 + 0) = 
isin(0/2)co(r/2 — 0) up to small corrections due to tunnelling between the wells. 
With an exact it phase shift, co(t) = for t > t/2 as in the non-interacting case. But 
for O 7^ 7r, both modes are nonlinearly coupled. Solving Eqs.(T5]) numerically, we obtain 
the results of Figf2](right) for the final population pq. Comparing to the full solution of 
the GPE, both curves show almost the same behaviour in the close vicinity of 6 = n, 
and agree qualitatively for other phase shifts: they both present a maximum and a 
minimum. This differs from the results of [16] . probably because there the eigenstates 
of the linear Schrodinger equation were used for the expansion. Moreover, a more 
accurate agreement could not be expected because the two-mode approximation is 
valid only when the low energy states in the separated wells are weakly modified by the 
many-body interactions. Estimates available in the literature, e.g. Ref.JB], show that 
the two- mode approximation is limited to an interaction strength g -C 2 3 / 4 y / 7r w 2.98, 
given our trap potential. 

In FigOTright), we compare the evolution of the weights |co,i(t)| 2 between the full 
numerical solution and the two-mode approximation for a phase shift close to tt. Again, 
the approximation reproduces the general features, with quantitative differences in the 
amplitudes and the frequencies of the oscillations. The curves also illustrate that for 
superpositions of instantaneous eigenstates, the adiabatic approximation fails due to 
the nonlinearity: the state after the phase imprint — it is almost identical to the lowest 
antisymmetric state — turns during the evolution into a state with a strong ground state 
admixture. This is due to the presence of an instability in the system. More precisely, 
the excitation spectrum of the antisymmetric state contains a soft mode [16j [23] . 
Solving the Bogoliubov-de Gennes equations, we have found that the first eigenmode 
acquires a purely imaginary frequency if the splitting d between the wells exceeds 
some critical value that becomes smaller with increasing interactions [jj The imaginary 
excitation frequency leads to an exponential growth of the ground state, and this 
instability is rooted in the Josephson effect [24, 25 : putting a relative phase different 
from 7r on the two arms of the interferometer, one introduces a tunnel current across 
the barrier that drives a population imbalance between the two wells. The interactions 
in each well subsequently enhance the phase difference, leading to a runaway effect for 

|| We were unable to reproduce quantitatively the imaginary eigenfrequency plotted in Ref. 16,. This 
may be related to the different ansatz for the perturbed order parameter used there. 
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the antisymmetric state. The instability is weaker at large splitting because tunnelling 
becomes exponentially suppressed. This is consistent with the behaviour of the growth 
rate found from the Bogoliubov analysis (see also [H]). For the chosen parameters 
the timescale of the process is actually dominated by the merging of the two wells 
that increases Josephson critical current. As illustrated in FigQJ right), the main 
population transfer occurs around the moment when the barrier between the two 
wells drops below the chemical potential. 

The previous results suggest that in a realistic setting, the condensate 
interferometer is very sensitive to fluctuations of the trapping potentials. It has also 
become clear recently that electromagnetic field fluctuations are particularly relevant 
for integrated atom optics in the neighbourhood of material microstructures held 
at a finite temperature (see ESI EZ] and references therein). We have simulated 
these fluctuations in terms of a random potential that fluctuates in time and space 
with spectral characteristics similar to those of thermal magnetic near fields: white 
noise and a lorentzian spatial correlation function EE]. We focus on the case that 
the correlation length is smaller than the distance between wells. This is realistic 
for microtraps whose splitting is larger than the distance to the underlying surface. 
Fig[3] shows the impact of increasing the noise strength for the interferometer output, 
averaged over a statistical sample of realizations. The steep structures around a 7r 
phase shift are smoothed, as expected. It is surprising that this happens already for 
very weak noise: if the scattering rate off the noise potential is denoted 7, significant 
changes are seen for jt ~ 0.07 already. The instability in the double-well potential 
thus makes the condensate more sensitive to thermal near fields. This differs from a 
non-split condensate in a single well, which is more robust in the presence of noise 
compared to an ideal g shown in 291. 

< Po > < Pi > < Pox > 




O/TT &/7Z &/7Z 



Figure 3. Average populations (po), (pi), and (p C x) (higher excitations) in 
a noisy interferometer for different values of the noise power spectral density: 
7 = 0.01 (dash-dotted line), 7 = 0.001 (dotted line), 7 = 0.0001 (dashed line), 
and without noise (solid line). The error bars indicate the statistical uncertainty 
of the simulation. 

Throughout this paper, we have used the mean field description based on the 
Gross-Pitaevski equation and neglected phase fluctuations. At fixed g, this can be 
justified for a sufficiently large number of atoms. The phase fluctuations in the initial 
condensate are still small even at finite temperatures < T < T v w N/(Q±/fl + 
(3(7/(4\/2)) 2 / 3 ) (in our units) [IT], using the regime of tight radial confinement. 
Similarly, the phase coherence of the split condensate 1F1 i 1 3| is maintained for operation 
times t < Tdiff w (2v / 3/g) 2 ^ 3 -/V 1 / 2 , following the procedure explained in Ref.[9]. 
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